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Abstract 

In this paper, we study in detail a perfect fluid cosmological model with time- varying "constants" 
using dimensional analysis and the symmetry method. We examine the case of variable "constants" 
in detail without considering the perfect fluid model as a limiting case of a model with a causal 
bulk viscous fluid as discussed in a recent paper. We obtain some new solutions through the Lie 
\ method and show that when matter creation is considered, these solutions are physically relevant. 
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I. INTRODUCTION 



Cosmological models with time- varying "constants" have been studied for quite some time 
ever since Dirac 1 proposed a theory with a time-varying gravitational constant G. Several 
works have investigated cosmological models with variable cosmological constant within a 
framework of dissipative thermodynamics 2-14 . In a recent paper 15 , we have studied a causal 
bulk viscous cosmological model with time-varying constants. We arrived to the conclusion 
that our cosmological model tends to a perfect fluid one in the matter predominance era. 
The approach of our work was to study the symmetries of the model beginning with the 
dimensional analysis of the field equations. The method of Lie group allowed us to arrive 
to the conclusion that under the hypotheses considered there is only one solution for our 
model, the scaling one that was trivially obtained using dimensional analysis. 

Since our viscous model tends to a perfect fluid one, the purpose of this work is to perform 
a detailed study of all the possible symmetries of a perfect fluid model with time varying 
constants showing that in this case it is possible to find more solutions in addition to the 
scaling one. In order to carry out this study, we begin in section 2 by outlining the equations 
that govern the model as well as the notation employed. In section 3, we review the scaling 
solution obtained in previous works highlighting the "assumed" hypotheses that we need 
to make in order to obtain a solution using dimensional analysis, these are: div(Tj) = 0, 
conservation principle, and that the relation G/c 2 remain constant for all value of t (cosmic 
time). We emphasize the special case uo = i.e. the dust case as this is the scenario that 
describes our model in the matter predominance era (in agreement with our previous paper). 
We discuss some interesting relationships that arise with the similarity method. 

In Section 4, we work towards finding other possible solutions to the field equations using 
the Lie group method. We start this section by rewriting the field equations in such a way 
that we can use the standard Lie procedure that allow us to find more symmetries. After 
outlining the equation and the constraint, we proceed to study some cases. The first one 
is the obtained previously using dimensional analysis since dimensional analysis is just a 
special class of symmetry. We would like to emphasize that the Lie method show us that 
one of the assumptions made with the dimensional method, G/c 2 = const., is at least correct 
from the mathematical point of view. This result allow us to validate completely the solution 
obtained through similarity. This solution connects perfectly with our previous work 15 
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In order to show that the other two solutions are physically relevant, we examine them 
in a matter creation scenario in section 5. We find that the horizon and entropy problems 
are solved when we consider matter creation and this leads to a set of physically relevant 
solutions. We note that no new constants or assumptions are required for this exercise. We 
conclude the paper by summarizing the results in section 6. 



II. THE MODEL 

We will use the field equations in the form: 

Rh - \gaR = ^j^ T ii + A (t)gij, (i) 

where the energy momentum tensor is: 

T ij = (P + P) u i u i - P9%j, (2) 

and p = up in such a way that u G [0, 1] . 

The cosmological equations are now (for a flat FRW Universe as the most recent obser- 
vations suggest us 16 " 19 .): 

2H> + 3tf 2 = -^p^ + c(i) 2 A(t), (3) 

3# 2 = ^p + c(t) 2 A(t), (4) 

where H = (/'/ /) is the Hubble function. Applying the covariance divergence to the second 
member of equation (1) we get: 

which simplifies to: 

pf + 3{u + l)pH = - P % + 4p-. (6) 

We assume that div(Tj) = which leads to the two following equations: 

p' + 3(oj + l)pH = 0, (7) 

A'c 4 G' . d 

p— + 4p- = 0. (8) 



ItxG G 
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Hence the field equations are: 



p' + 3(uj + l) P H 

A'c 4 a d 

h 4— 

87rpG G c 



2H' + 3H 2 



3H 2 




(10) 

(11) 
(12) 



(9) 



0. 



III. REVIEW OF THE SCALING SOLUTION 

As have been pointed out by Carr and Coley 20 , the existence of self- similar solutions 
(Barenblatt and Zeldovich 21 ) is related to conservation laws and to the invariance of the 
problem with respect to the group of similarity transformations of quantities with inde- 
pendent dimensions. This can be characterized within general relativity by the existence 
of a homotetic vector field and for this reason one must distinguish between geometrical 
and physical self-similarity. Geometrical similarity is a property of the spacetime metric, 
whereas physical similarity is a property of the matter fields (our case). In the case of per- 
fect fluid solutions admitting a homotetic vector, geometrical self-similarity implies physical 
self-similarity. 

As we show in this section as well as in previous works, the assumption of self-similarity 
reduces the mathematical complexity of the governing differential equations. This makes 
such solutions easier to study mathematically. Indeed self-similarity in the broadest Lie 
sense refers to an invariance which allows such a reduction. 

Perfect fluid space-times admitting a homotetic vector within general relativity have been 
studied by Eardley 22 . In such space-times, all physical transformations occur according 
to their respective dimensions, in such a way that geometric and physical self-similarity 
coincide. It is said that these space-times admit a transitive similarity group and space- 
times admitting a non-trivial similarity group are called self-similar. Our model i.e. a flat 
FRW model with a perfect fluid stress-energy tensor has this property and as already have 
been pointed out by Wainwright 23 , this model has a power law solution. 

Under the action of a similarity group, each physical quantity <fi transforms according to 
its dimension q under the scale transformation. For space-times with a transitive similarity 
group, dimensionless quantities are therefore spacetime constants. This implies that the 
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ratio of the pressure of the energy density is constant so that the only possible equation of 
state is the usual one in cosmology i.e. p = up, where u is a constant. In the same way, the 
existence of homotetic vector implies the existence of conserved quantities. 

In this section we would like to review the solution obtained through Dimensional 
Analysis 24-25 . Our starting point is the condition div (T/) = which allows us to obtain 
one of the dimensional constants that we need in order to apply the method of dimensional 
analysis. Therefore, from eq. (11) 

p' + 3(p + p)H = 0, (13) 

we obtain the following relation between the energy density and the scale factor as well as 
the constant of integration that we shall need for our subsequent calculations: 

p = AJ~'^ +l \ (14) 

where is the integration constant that depends on the equation of state that we want to 
consider i.e. constant u, \AJ\ = L 3 ^ 1 ^ 1 MT~ 2 . We consider a second dimensional constant 
by considering the relation G/c 2 = B, where, B is the constant. This is a hypothesis which 
is necessary in order to apply dimensional analysis. In this next section, we will see that this 
condition on G and c is mathematically correct. Our purpose here is to show that no more 
hypotheses are necessary to solve the differential equations that govern the model (see 26 ~ 28 
for standard text-book on Dimensional Analysis and 29 for applications to Cosmology). 

Therefore, if we take into account the standard dimensional procedure, we find that the 
set of governing parameters are M = M {A^, B, t} , which bring us to obtain the next 
relations: 

_2_ 2 . i 2(1 — A 

G oc AZ +1 B—^ + H^- , 

. - TT 1 (1-7) 

c oc AZ + B^t^+i~, 
p oc B'H- 2 , 

f oc A^BTfit^, (15) 
k B 6 oc A^B^" 1 ^, 

-^2- -2 -4 

A OC AZ +1 Br&triZ , 

7-1 

where 7 = 3(u + 1) — 1, and q is the deceleration parameter. 

From the set of equations (15), we see that ^ = B (trivially), / = ct (the horizon problem 
is missing in this model), A oc -^2 oc f~ 2 for all value of 7 i.e. Vcj. 
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We would like to point out that the same results have been obtained by Midy and Pettit 
in a more general context 30 . 

As we have indicated in the introduction, in a recent paper 15 we have studied a full causal 
bulk viscous cosmological model with time varying constants. The main conclusions in that 
paper was that under the "assumed" hypotheses, i.e. div(Tj) = 0, conservation principle, 
and IT oc p, the bulk viscous pressure behaves as the energy density and we obtained the 
following results: 

1. there is only one solution for the field equations (scaling solution), 

2. the bulk viscous pressure behaves as an adiabatic matter mechanism (solving the 
entropy problem), 

3. the "constants" G, c and A are decreasing functions on time (solving the horizon 
problem) and 

4. in the matter predominance era, i.e. uj = 0, we cannot consider the bulk viscosity and 
therefore our viscous fluid tends to a perfect fluid one (see 15 for details). 

Due to these conclusions, in the present paper, we have studied a perfect fluid model with 
variable "constant" emphasizing the case of matter predominance era or a "dust" solution; 
in other words, considering u = in the equation of state. In this case, we obtain the 
following solutions: 

G oc t' 2/3 ,c oc t~ 1/3 ,p oc t~ 2 

f oc t 2/3 ,#oct°,Aocr 4/3 ,g= 1/2 (16) 

while in the viscous model, the solutions are: 

G oc AS^k^r 4 -^), (17) 

c oc A^k^r 1 ' 1 ^, (18) 

p oc k- b ~ 1 t b ~ 1 oc IT, (19) 

/ oc A^kf^r^), (20) 

A oc i^fcf^t^, (21) 
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where, a = 3(u>+ 1 + x) — 1, b = 7 — l,(see 15 for details). If we consider 7 = 1/2, uj = and 
x = (viscous pressure vanishes and therefore k y vanishes too) then we obtain the same 
results as in equation (16), but, as we have indicated previously this approach is inconsistent 
when u = and for this reason we have studied the perfect fluid case. 

As we can see in the perfect fluid case, "constants" G, c and A are decreasing functions 
on time, but in this case decrease slower than in the radiation predominance era, while p 
and / behave as in the FRW model. 

To obtain these solutions we needed two dimensional constants, viz., and B. In the 
next section, we shall show that constant B is a reasonable hypothesis since with the Lie 
group method such condition holds as a result in the scaling solution. 



IV. LIE METHOD 



As we have seen earlier, the n — monomia is the main object in dimensional analysis. It 
may be defined as a product of quantities which are invariant under changes of fundamental 
units. 7r — monomia are dimensionless quantities, their dimensions are equal to unity. 
Dimensional analysis has the structure of a Lie group 31 . The it — monomia are invariant 
under the action of the similarity group. On the other hand, we must mention that the 
similarity group is only a special class of the mother group of all symmetries that can be 
obtained using the Lie method. For this reason, when one uses dimensional analysis, only 
one of the possible solutions to the problem is obtained. 

As we have been able to find a solution through dimensional analysis, it is possible that 
there are other symmetries of the model, since dimensional analysis is a reminiscent of 
scaling symmetries, which obviously are not the most general form of symmetries. Hence, 
we shall study the model through the method of Lie group symmetries, showing that under 
the assumed hypotheses there are other solutions of the field equations. In this section we 
shall show how the lie method allows us to obtain different solutions for the field equations. 
In particular we seek the forms of G and c for which our field equations admit symmetries 
i.e. are integrable (see 32_37 ). 

An alternative use of the Lie groups have been performed by M. Szydlowki et. al. 38-39 
where they study the Friedman equations in order to find the correct equation of state 
following pionerr works of Collins 40 . 
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In order to use the Lie method, we rewrite the field equations as follows. From (9) 
(10), we obtain 



2 

'J 



= "— (P + P)> (22) 



and therefore 



2W = — (p + p). (23) 



c 2 



From equation (11), we can obtain 



# = „ xx , (24) 

3 w+lp' 1 } 



therefore 



127i(u + iy -p. (25) 
P J c 



Taking 12n (u + l) 2 = A and then expanding, we obtain 



n P' 2 , A G 



P"= r - + A- 2 p\ (26) 
p c z 

Now, we apply the standard Lie procedure to this equation. A vector field X 

X = at,p)d t + r){t,p)d p , (27) 

is a symmetry of (26) iff 

-e/t - vf P + % + (2^ P - e«) p' + (^ - 26 P ) p /2 - w 3 + ... 

... + ( Vp - 26 - 3p'£ P ) / - fr t + (tj p - 6) p' - P%] fp> = 0. (28) 

By expanding and separating (28) with respect to powers of p', we obtain the overdetermined 
system: 

U + P _1 £p = 0, (29) 
Vpp -2i tp + p- 2 r l -p~ l r lp = 0, (30) 
2vt P ~ & - 3A^p 2 ^ p - 2p~ l r lt = 0, (31) 

% -a(%- 2G^) p 2 i - 2 V A%p + (rj p - 26) A^p 2 = 0. (32) 
Solving (29-32), we find that 

p) = -2et + a, ri(t, p) = (bt + d) p, (33) 



2- => — = 5 = const. (35) 



subject to the constrain 

with a, b, e, and d as constants. In order to solve (34), we consider the following cases. 



A. Case I: b = and d - 4e = 

In this case, the solution (34) reduces to 

G' _ d G 
~G~ 2 ~c^d 

which means that "constants" G and c vary but in such a way that the relation § remains 
constant. 

The solution obtained through Dimensional Analysis needs to make this relations as 
hypothesis in order to obtain a complete solution for the field equations. This case shows 
us that such hypothesis is correct (at least has mathematical sense). 

The knowledge of one symmetry X might suggest the form of a particular solution as an 
invariant of the operator X i.e. the solution of 

dt dp 



v(t, P y (36) 

this particular solution is known as an invariant solution (generalization of similarity solu- 
tion), therefore the energy density is obtained as 

P = — 3» ( 37 ) 



-2et + a Aep r (2et - a) 
for simplicity we adopt 

P = Pot' 2 , (38) 
Once we have obtained p, we can obtain / (the scale factor) from 

p = 4,r 3(w+1) =>f = (A,t)WW , (39) 

in this way we find H and from eq. (10), we obtain the behaviour of A as: 

c 2 A = 3H 2 - ^fp, (40) 

and therefore, 

A = (3/? 2 - 8nB Po ) j^ = ^ 2 . (41) 



If we replace all these results into eq. (12), then we obtain the exact behaviour for c, i.e., 
where A 



SnBpo , with A G M + , i.e. is a positive real number and thus, 

c = c r a , 

with a = (^) . 

Hence, in this case we have found that (see fig.l): 

G = G t~ 2a , c = c t~ a , A = A r 2(1 - a \f = (AJ)^ , p = p t- 2 . 



(43) 



(44) 



This is the solution that we have obtained with dimensional analysis in the previous section 
and we shall show this is only solution compatible with our previous solution 15 obtained in 
the framework of a bulk viscous fluid (full causal theory). 





FIG. 1: We see the behaviour of G, c and A for the first class of solutions for different values of 
a : a = 0.5 (solid curve), a = 0.9 (dotted curve), a = 0.3 (matter era)(dashed curve). In all cases 
the constants are decreasing functions. 



B. Case II, b = a = 



In this case, we find that 

% = St", (45) 



c 2 



where x = S — 2 and 6 = O n following the same procedure as above, we find that 

£ = Po t-\ (46) 
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we must impose the condition signed) = sign(e), i.e., 6 E IR + , in order that the solution 
has some physical meaning that the energy density is a decreasing function of time t. It is 
observed that if d — 4e then we obtain same solution that the obtained one in the case I. 
The scale factor is found to be 

f = K f t^+ T >, (47) 
where Kf is an integration constant, and therefore, the Hubble parameter is: 

H = , - x , (48) 

which is similar to the scale factor obtained in case I. To obtain the behaviour of the 
"constants" G, c and A, we follow the same steps as in case I, i.e., from 

c 2 A = 3tf 2 - *H£.p, (49) 



c 2 



we obtain the behaviour of A being: 

A = ^ (5°) 
where, / = {K\ — K 2 ), K\ = -. ^ , 2 and K 2 = 8irpoB i.e., I E 1R + . Therefore, 



If we substitute all this results into the next equation 



A'c 4 G' d 

8^ + G- 4 7 = °' (52 » 

we obtain an ODE for c, i.e., 

- (A - 2) = - (A - 2 + 8) - (53) 
c i 

where, A = f 1 —~ , A S R~, which leads to 

c = c t~ a (54) 
with a = (l + T^g) scuh that a E [0, 1). In this way we can find the rest of quantities: 

G = G t- 2(Q+1)+<5 , A = A r 2(1 ~ a) , (55) 

note that a < 1. The case a = 1 <^=^ 5 = is forbiden and a = brings us to the limiting 
case of the G, A variable cosmologies 41 . 

We notice that this solution is very similar to the case I but in this case all the parameters 
are perturbed by 5 and more important is the result, ^ = Bt* (see figs. 2, 3 and 4). 
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0.5 



FIG. 2: Time variation of c(t) for the second class of solutions for different values of a : a = 0.5 
(solid curve), a = 0.9 (dotted curve), a = 0.1 (dashed curve) and a = 0.000001 (long dashed 
curve), the last solution describes the case c(t) = const. 

C. Case III, b = e = 



Following the same procedure as above, we find in this case that such restrictions imply 
£(t, p) = a, r)(t, p) = dp and therefore: 



which brings us to: 



9L-2--- 

G c a 



% = Kexp(-at), 
c z 



where ^ = a and note that [K] = [B] i.e has the same dimensional equation, 

dt dp dt dp 

4(t,p) v{t,P) a dp 

this expression only has sense if a G IR~, note that [a] = T _1 . 

The scale factor / satisfies the relationship: 



p = 4,/" 3(w+1) =>f = K, exp (at) TO 



(56) 



(57) 



(58) 



that is to say, it is a growing function without singularity In this way, we find that 

a 



H 



3(u + l] 



conts. H > 0. 



(59) 



(60) 
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t 

FIG. 3: The variation of the gravitational "constant" G(t), for different values of a and 5. : a = 0.5 
and 5 = 1 (solid curve), a = 0.9 and 5=1 (dotted curve), a = 0.1 and 5 = 1 (dashed curve) and 
a = 0.000001 and 5 = 5 (long dashed curve), the last curve describes a growing solution. 



The cosmological "constant" is obtained as 



a 2 

c 2 A = 2 - 8nKp =^ c 2 A = I, (61) 

3 (u + 1) 

note that [/] = T -2 , if we replace all these results into eq. (12) then we shall obtain the 
exact behaviour for c, i.e. 

+ 2 ) - = a, (62) 



J&ixKpQ J c 
and hence, 

c = K exp(cot) , (63) 
where Cq = -, — p r- with cq G IR~ since a G 1R~, that is, c is a decreasing function on time 

t. 



In this case, we have found 



c = i^exp(c t), (64) 
G = G exp((-a + 2c )t), (65) 
A = Zexp(c t)" 2 , (66) 



13 




- 



0.5 1 1.5 2 

t 

FIG. 4: Time variation of A(t) for the second class of solutions for different values of a : a = 0.5 
(solid curve), a = 0.9 (dotted curve), a = 0.1 (dashed curve) and a = 0.000001 (long dashed 
curve). In all cases, A(t) is a decreasing function. 

therefore the solutions for this case are (see fig. 3): 

G = Go exp((— a + 2co) t), c = K exp(cot), A = / exp(cot) -2 , (67) 
p = po exp (at) , / = K f exp (at)^ry . (68) 




FIG. 5: We see the behaviour of constants for the third class of solutions for different values of 
Co ■ Co = —0.5 (solid curve), Co = —0.3 (dotted curve), cq = —0.01 (dashed curve) and a = —0.1. 
In all cases, A(t) is a growing function. 
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As we can see, equation (34) allows us to obtain more cases, however, we examine their 
physical validity by considering matter creation in our model in the next section. The 
behavior of this class of solutions for different values of Co is shown in figure 5. 



V. MATTER CREATION 

As we have seen in the previous section, solutions IV B and IV C suggest us a new 
scenario since these solutions may describe the early universe in a very different way than 
our previous solution 15 . For this reason, in this section, we study briefly the important case 
in which adiabatic matter creation 42-45 can be taken into account, in order to get rid of the 
entropy problem. The matter creation theory is based on an interpretation of the matter 
energy-stress tensor in open thermodynamic systems, which leads to the modification of the 
adiabatic energy conservation law and as a result including the irreversible matter creation. 
The matter creation corresponds to an irreversible energy flow from the gravitational field to 
the constituents of the particles created and this involves the addition of a creation pressure 
p c in the matter energy-momentum tensor which we discuss below. 

The field equations that now govern our model are as follows: 

2H' + 3H 2 = -^^(p + Pc ) + c 2 (t)A(t), (69) 

3R2 = ^f P + c2(t)A(t) ' (70) 
n' + 3nH = if), (71) 

and taking into account our general assumption i.e. 

, (4c d G A , c^A, 
^ V c G ) * + 8nG _U ' {U) 

with T-.j = 0, we obtain the two equations 

p' + 3(p + p + p c )H = 0, (73) 

and 

A ' c4 G ' A C ' n 

— — + — - 4- = 0, 74 

87rGp G c 

where n is the particle number density, ip is the function that measures the matter creation, 
H = f'/f represents the Hubble parameter (/ is the scale factor that appears in the metric), 
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p is the thermostatic pressure, p is energy density and p c is the pressure that generates the 
matter creation. 

The creation pressure p c depends on the function ip. For adiabatic matter creation this 
pressure takes the following form: 



Pc 



3nH 

The state equation that we next use is the well-known expression 



(75) 



p = up, (76) 

where u = const, and u 6 (— 1, 1]. We assume that the matter creation function follows the 
law 45 : 

ip = 3f3nH, (77) 

where (3 is a dimensionless constant (if f3 — then there is no matter creation since — 0). 
The generalized principle of conservation T- j = 0, for the stress-energy tensor (73) leads us 
to: 

f/ + 3(u + l)(l-P)pH=0. (78) 
Therefore, the new set of field equations are: 

2H' + 3H 2 = - 8 -^(p + Pc ) + c 2 (t)A(t), (79) 
c\t) 

p' + 3(cu + 1) (1 - p) pH = 0, (81) 

a'c 4 a d 

8^ + G- 4 7 = °- (82) 
Now, using the same procedure that was followed in section 4, we find the new equation 
on which we can apply the Lie method. We rewrite the field equations as follows: from (79) 
— (80) we obtain 

and therefore 



2(H)' = -^f(p + p c + P ). (84) 



From (81), we can obtain 



H --3{u + i){i-py (85) 



1(3 



therefore 




(86) 



making 12-n (u + 1) (1 — f3) 2 = A and expanding, we obtain 




(87) 



Therefore, we obtain the same equation that we obtained in section 4, namely equation 
(26) except the constant A. Hence, within this framework the entropy and horizon problems 
are solved for our model. 

VI. CONCLUSIONS. 

In this paper we have studied the behaviours of time- varying "constants" G, c and A in a 
perfect fluid model. We began reviewing the scaling solution obtained through dimensional 
analysis. He have shown that this solution connects with our previous solution 15 where we 
studied the behaviour of the "constants" G, c and A in a full causal bulk viscous model 
arriving to the conclusion that this model tends to a perfect fluid one when we impose the 
condition uj = in the equation of state (dust solution). 

To obtain this solution, we imposed the assumption, div(Tj) = 0, from which we obtained 
the dimensional constant A w that relates p oc f- 3 ( u}+1 ) and the relationship G/c 2 = const. = 
B remaining constants for all value of t, i.e. G and c vary but in such a way that G/c 2 
remain constant. With these two hypothesis, we have obtained the scaling solution that 
connects perfectly with one obtained in our previous paper 15 , i.e. with the bulk viscous 
model in the matter-dominated era. 

In this context, the solution obtained through dimensional analysis show us that the 
"constants" G, c and A are decreasing functions of time, but in this case decrease slowly than 
in the radiation predominance era, while p and / behave as in the FRW model solvingthe 
horizon problem. 

Since we have been able to found a solution through similarity, i.e. through dimensional 
analysis, it is possible that there are other symmetries of the model, since dimensional 
analysis is a reminiscent of scaling symmetries, which obviously are not the most general 
form of symmetries. Therefore, we studied the model through the method of Lie group 
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symmetries, showing that under the assumed hypotheses, there are other solutions of the 
field equations. 

The first solution obtained is the already obtained one through similarity, but, in this case 
we have showed the condition G/c 2 arises as a result and not as an ad-hoc condition. We 
also have studied two other cases which can be considered as physically relevant solutions 
since / is a growing function on time and p is a decreasing function on time. They could 
describe very early cosmological solutions (inflationary ones) but in this context we cannot 
solve the entropy problem. We have considered matter creation and and using a recasted 
set of equations, shown that when matter creation is taken into account, the horizon and 
entropy problems are solved for the two solutions. 
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